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Abstract 

a ■ 

The electronic structure of interfaces between lattice-mismatched semicon- 
C I ductor is sensitive to the strain. We compare two approaches for calculating 

such inhomogeneous strain - continuum elasticity (CE, treated as a finite 
difference problem) and atomistic elasticity (AE). While for small strain the 
£Sj ' two methods must agree, for the large strains that exist between lattice- 

mismatched III-V semiconductors (e.g., 7 % for InAs/GaAs outside the lin- 
' earity regime of CE) there are discrepancies. We compare the strain profile 

obtained by both approaches (including the approximation of the correct C2 
symmetry by the C4 symmetry in the CE method), when applied to C2- 

^-j ■ 

q ' symmetric InAs pyramidal dots capped by GaAs. 
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I. INTRODUCTION 



One of the leading methods for growing semiconductor quantum dots is via the controlled 
coarsening of a film of a material that is strained with respect to the substrate on which it 
is grown. EH This ("self-assembled") coarsening/roughening is a result of lattice-mismatch- 
induced strains. The dots are often capped by the substrate material, thus extending the 
strain around the dot to all angular directions. Not surprisingly, the interpretation of the 
electronic structure of such dots is profoundly affected by their strain profile. Thus, in 
order to calculate or interpret the measured electronic structure, one has first to calculate 
or measure the position dependent strain tensor e a p. 

The three basic approaches to calculating such strains are: 

(i) Harmonic Continuum Elasticity: Here, one uses classical elasticity! within the har- 
monic approximation. For a cubic system, the strain energy per atom, Eqe, is 



Ece - -^Cu(t 2 xx + t 2 yy + t 2 zz ) + — C u (e 2 yz + e 2 zx + e 2 xy ) + VC 12 (e yy e zz + t zz e xx + e xx e yy ), (1) 



where V is the equilibrium volume, Cy are cubic elastic constants and e a p is the strain 
tensor. We illustrate the predictions of harmonic continuum elasticity for a 2D film, since 
this is going to be used as a test case. In the absence of shear strain (e a p oc 5 a p), for a fihn 
coherently grown on a substrate with parallel lattice constant a s , the strain components are 



£|| — £xx — Cyy 



^eq 



e ± = e zz = (2) 

Q>eq 

where a eq is the equilibrium lattice constant of the unstrained material and c is the per- 
pendicular lattice constant of the strained film. The equilibrium value of this c— axis is 
determined from dE C E/de± = 0, yielding, 

CegKG) -l = [2-3g(G)] (3) 

^eq 

where the "epitaxial strain reduction factor" for orientation G of the c— axis is 



and A = C44 — 1/2 (Cn — C12) is the elastic anisotropy, £> = 2/3 (Cn + 2Ciq) is the bulk 
modulus and 7(G) is a purely geometric factor given in Ref. p|. For principal directions, 
7(001) = 0,7(011) = 1 and 7(111) = 4/3. Equations (|)-(D are use d routinely to predict 
tetragonal distortions of strained films!. The harmonic continuum elasticity method has 
been recently applied to pyramidal quantum dots by Grundmann et aM and by Pryor et 
al® 

(ii) Atomistic elasticity: Here, one avoids a continuum description and describes the 
strain energy in terms of few-body potentials between actual atoms 

Eae = E ^(R< - a,-) + E v 3 (e ijk ) + • ■ ■ , (5) 

ij ijk 

where V2 is a two-body term, V3 is a three-body function of the bond angle, 0y&. The 
functional form of these terms is taken to be strain-independent. The strain is determined 
by minimizing Eae with respect to atomic positions {R}. Like the continuum elasticity 
approach, only the cubic elastic constants are used as inputffii. However, unlike the CE 
approach, here (a) optical phonon modes can be described,ifl (b) harmonicity is not assumed, 
and (c) the atomic level symmetry is retained. The last point is illustrated in Fig. 1 that 
shows a regular pyramid representing the quantum dots from experiments of Grundmann 
et a/.i In a continuum representation, the strain is equal on the {110} and {110} facets, 
while in an atomistic description these two facets can have different strain if the pyramid is 
made of a zincblend material. Atomistic elasticity has been widely used to determine strain 
in alloys,i0 superlattices0 and dots^'tl where V2 and V3 of the Eq. (5) are taken from 
Keating's valence force fieldSSi (VFF) model. 

(Hi) Atomistic quantum-mechanical approach: Here one does not have to assume any 
model for inter-atomic interactions as in the atomic elasticity. Instead, one explicitly com- 
putes the total electron and nuclear energy £7 tot [{Rj}] for each atomic configuration {R;} 
directly from a quantum-mechanical Schrodinger equation. Atomic symmetry is retained 



and harmonicity is not assumed. This approach has been used for small (< 100 atom) wires 
and clustersEHHI but it impractical for ~ 100 A dots (~ 10 5 atoms). 

The three approaches to the calculation of strain - harmonic continuum elasticity, (an- 
harmonic) atomistic elasticity, and the atomistic quantum mechanical approached - have 
been recently compared for InAs/GaAs strained superlattices.0 However, no comparison 
exists for 0-dimensional quantum dots. Here we perform parallel calculations for the strain 
e(R) of a pyramidal InAs dots (Fig. 1) surrounded by GaAs using the two approaches that 
are practical for large dots: continuum elasticity and atomistic elasticity. We find that: (i) 
the strain profiles obtained via continuum elasticity are in qualitative agreements with those 
found by atomistic elasticity; (ii) the atomistic elasticity produces different strains on the 
two facets ({110} or {110}) of the Zincblend pyramidal dots (see Fig. 1), corresponding to 
the physical C2 symmetry, while the continuum elasticity approximates this as C4 symmetric 
strain, (iii) the quantitative discrepancy resides mostly inside the dots, while the difference 
in the barrier region is smaller; These differences are traced back to the fact that the strain 
lies outside the domain of validity for the linear elasticity. We illustrate this point by con- 
trasting the predicted c eq (a s , G)/a eq ratio of coherent 2D films, as obtained by harmonic 
continuum elasticity [Eq. (3)] and atomistic elasticity. Differences are noticeable already for 
1 % biaxial strain, whereas the controlled-coarsening ( "self-assembled" ) growth method for 
quantum dots needs to deal with larger mismatches (7 % for InAs/GaAs and InP/GaP). 
Finally, the consequences on the electronic structure of the different strain profiles obtained 
for dots using CE and AE are illustrated. 



II. METHODS OF CALCULATIONS 

A. Continuum elasticity for dots 

In the CE approximation the strain is determined by minimizing the elastic energy given 
in Eq. 1. To account for the lattice mismatch we assume the coordinates are fixed to the 
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barrier material, and treat the island as expanded barrier material (with different elastic 
constants). This is accomplished by the modification 

E CE -> E CE - a(r)(e xx + e yy + e zz ) (6) 

!0 barrier 
(7) 
(Cxxxx + 2C xxyy ) (a/ -a B )/a B island 

where a B and aj are the unstrained lattice constants for the barrier and island material 

respectively, and C xxxx and C xxyy are the elastic constants for the island material. A piece 

of island material with no external forces acting on it will have its energy minimum shifted 

to 

6ij = 5ij{ai - a B )/a B . (8) 

This fictitious strain corresponds to unstrained island material and must be subtracted. The 
corrected strain is still computed with derivatives in the barrier's coordinates and must be 
converted to the the island coordinates through multiplication by Thus, the physical 



strain is given by 



e phys = ^£( e .. _ g.^ ai _ a b )/a B ) (9) 
ai 

where = + ^r)/2 is the strain computed directly from the displacement Ui which 
minimizes Eqe- 

A numerical solution requires some kind of discretization. We define the displacements 
Ui on a cubic grid, thereby maintaining the cubic symmetry of the crystal. The strain is 
expressed in terms of forwards or backwards differences by 

4 = (Afuj + A±«j)/2 (10) 

I n * I 

A-/(r) = /W - (12 ) 
where is the lattice vector in the i direction. Symmetric differences 
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A|/(r) 



/(r + n t )-/(r-n,) 

2|n. 4 | 



(13) 



(14) 



are undesirable since they give unphysical low energy configurations which oscillate with 
period 2 1 | . For example, a displacement u x {r) = sm(iix/\h x \) has e xx = when constructed 
using symmetric differences. The oscillatory solutions cannot be simply discarded since 
they mix with the physical ones. Non-symmetric derivatives are also problematic since a 
particular choice will single out a direction in space. The solution is to average Eqe over 
all permutation of ± on each of the three difference operators. That is, we take (E +++ + 
E~ ++ + E + ~ + + ...)/8. Physically this corresponds to taking the energy density at each site 
to be the average of the energy densities from each adjoining octant. 

The elastic energy is a quadratic function of the displacements, which is easily minimized 
using the conjugate gradient algorithm. For the barrier material the strains are computed 
directly using differences (now there is no impediment to using symmetric differences). In 
the island material we then apply the correction in Eq. (^|). 



In the VFF model, the strain energy is expressed as a functional of atomic positions, 
{Ri}, as 



Here, d?- denotes the ideal bond length between atoms i and j, and 8 is the ideal bond angle. 
For the Zincblend structure, cos#o = —1/3. The local-environment-dependent coefficients, 
aij and Pijk, are fitted to the elastic constants of bulk materials^. The long-range Coulomb 



B. Atomistic valence force field for dots 



Ejus = E ^(R* - R*) + £ V 3 {9 ijk ) 



ij ijk 




[(Rj - Ri) ■ (R k - Ri) - cos0 o <-4] 2 - 



(15) 
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interactions of Ref. p\ are neglected which causes a slight deviation from the measured bulk 
properties!. In this case, the elastic constants of a pure bulk Zincblend material are given 



as 



4r 

G\\ — C\2 = — (3 

4r a + p 

where r is interatomic bond length. Because Eq. ( P^|) contains only two free parameters, it 
is impossible to fit three arbitrary elastic constants. Nonetheless, for zincblend materials a 
and (3 may be chosen so the C"s fit within a few percent of the measured values. Table I 
gives the elastic constants of bulk GaAs and InAs calculated from Eq. (|l^) using en's and /3's 
of Ref. P]. The elastic constants obtained differ a bit from the experimental values (since 
Coulomb corrections to Eq. ( |i~6|) are neglected), but we will use them consistently in both 
our continuum elasticity and atomistic elasticity studies. For purposes of comparison, Table 
I also contains the experimental elastic parameters. Note that the VFF method with the 
standard parameterization of Eq.(16) reproduces well the values and trends in the formation 
enthalpies of strained GaP/InP structures, as obtained from first-principlesi. 

The relaxed atomic configuration is obtained by conjugate gradient minimization^ of 
Eae with respect to the atomic positions. At each minimization step, the atoms are displaced 
along the conjugate direction {h} by a finite increment A, as Rj —>■ Rj + Ahj. A line 
minimization of Eae along the conjugate gradient direction to find A that minimizes Eae 
is done by taking advantage of the fact that Eae is a fourth-order polynomial that depends 
on only the relative positions, Rj — Rj, of each atom: 

E AE [{R + Ah}] = EAE[{Ri - Rj}] 

+A£«[{R 4 - R,}, {h, - h,}] + A^KR, - R,}, {h, - h,}] 
+A 3 £( 3 )[{R - R,}, {h, - h,-}] + A 4 £( 4 )[{R, - R,}, {h, - h,-}]. (17) 
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The energy terms, E^-\E^ 2 \E^> and E^\ are also simple polynomials of {R} and {h}. 
The increment A minimizing the elastic energy is then obtained by solving exactly 

0^[{R+Ah}] = E(1) + 2XE(2) + 3A2jE(3) + 4A 3 E4 = 

OA 

Figure 2 illustrates how the local strain is calculated. After the atomic positions are re- 
laxed by minimizing Eae, the local strain tensor e at a cation site (for cation- mixed systems) 
is calculated by considering a tetrahedron formed by four nearest neighboring anions. The 
distorted tetrahedron edges, Ri 2 , R23 and R 34 , are related to the ideal tetrahedron edges 
R? 2 ,R^3 and K° u via 



Rl2,x 


R23,x 


R3A,x 


( 


1 ^xx 


tyx 


^zx 


\ 


ri '12,a; 


f?0 

^23^ 


-"-34,0; 


R\2,y 


Rz3,y 


Rs4,y 




£xy 


1 e yy 






R° 


R° 

23,y 


R° 

34,1/ 


y Rl2,z 


R23,z 


-^34,2 j 


V 


^xz 


e yz 


1 + e zz 


J 


y -"-12,2 


-^23,2 


-"-34,2 J 



(19) 



The ideal tetrahedron edges are {R } = {[110]a/2, [011]a/2, [I10]a/2}, where a denotes the 
equilibrium lattice constant of the cation, i.e., acaAs for Ga atoms and ai n A s for In atoms. 
The local strain, e, is then calculated by a matrix inversion as 



^xx 


£yx 


^zx 


\ 




Rl2,x 


R23,x 


R34,x 


\ 


( f?0 

-"-12,0; 


Ix 23,x 


-"-34,2; 


\ 


e xy 


e yy 


e zy 






Rl2,y 


R%3,y 


R%4:,y 






""-23,2/ 


-"-34,1/ 




^ ^xz 


£yz 


^zz 


J 




y Rl2,z 


R23,z 


-^34,2 


J 


dO 

y -"-12,2 


DO 

-"-23,2 


RO 

-"-34,2 


/ 



- I, (20) 



where / is the unit matrix. 



III. RESULTS 

A. Comparison of strain profiles 

Figure 3 shows e xx ,e zz and Tr(e) = e xx + e yy + e zz as obtained by continuum elasticity 
(dashed lines) and by atomistic elasticity (solid lines) as a function of the position from the 
pyramidal center along the [110] direction at a height z = h/3 from the base (see Fig. 1). 
The corresponding differences in strains, Ae = e{CE) — e(AE), are given as the solid lines 
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in Figure 4. We note that the grid points of the continuum elasticity calculation are chosen 
to be commensurate with the cation positions of the ideal GaAs Zincblend structure for 
consistent comparisons of the two approaches. The largest differences occur around the 
interfaces between the dot and the cap. A significant discrepancy is also found inside the 
quantum dot where the InAs experience large compressive strains: e xx of the continuum 
elasticity is found to be more compressive than that of the atomistic elasticity, while the e zz 
of the CE is more tensile. A similar comparison is given in Fig. 5, but this time the position 
vector is along the Z = [001] direction, starting from the substrate, going through the 
wetting layer into the pyramidal tip and then into the capping layer. Again, the discrepancy 
is largest around the interfaces, while the strains in the barrier (GaAs substrate and capping 
layer) agree within 0.5%. 

Figure 4 illustrates the extent to which the continuum elasticity description misses the 
correct atomic symmetry. In a pyramid made of Zincblend materials on the (001) substrates, 
the {110} and {110} facets are symmetrically inequivalent (Fig. 1). Indeed, the atom- 
istic calculation produces different strains. The dashed lines in Fig. 4 show the difference 
e^ E ([110]) — e^QllO]) for these two directions. We see that the anisotropy is pronounced at 
the interfaces. For the atomistic elasticity calculation, we construct the pyramidal structure 
to have an In atom at the pyramidal tip. This tip In atom has: (i) two-As atoms that 
belong to the InAs dot and lie along the [110] direction; and (ii) the other two-As atoms 
that belong to the GaAs capping layer and lie along the [110] direction. Considering only 
the local strain of the tip atom, one expects larger compressive e xx and e yy along the [110] 
direction than along the [110] direction, based on the atomic configuration. By the same 
token, the atoms at the {110} interfaces experience larger compressive e xx and e yy , than 
those at the {110} interfaces, for this particular choice of the pyramidal geometry. In the 
continuum description, e^ E ([110]) = e^QllO]) and this effect is missing. 
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B. The origin of the differences — a simple test case 



We know that the continuum and atomistic models, starting from the same input elastic 
constants, must agree in the limit of small strain and a large system. To study the rate 
at which the two methods diverge with increasing strain, we consider the simple case of 
biaxial strain. As Eqs. (FJ)-(^) show, for a 2-dimensional film that is constrained on a (001) 
substrate, continuum elasticity predicts 

2C 12 



en 



C Q>eq 
(X s Cl e q 



CE 



(21) 



Figure 6(a) compares this result with that obtained via atomistic elasticity, as a function of 
the relative film/substrate mismatch e\\ = (a s — a eq )/a eq . Similarly, for the (110) strain, 

_ Cn + 3C*i2 — 2C44 ^22^ 

£|| Cn + C\2 + 2C44 

and the corresponding comparison of the continuum and atomistic elasticity is shown in 
Fig. 6(b). We see that the discrepancy rises linearly, reaching 4 % for a lattice mismatch 
of 7 %, characteristic of InAs/GaAs. This difference is comparable to that found between 
CE and AE around the interfaces of the quantum dots (Figs. 3-5). Thus, the discrepancy 
simply reflects the departure from the linearity regime of the continuum elasticity. 



C. Consequences of the different strains in continuum elasticity and atomistic 

elasticity 

The existence of different strain magnitudes and even symmetries in a continuum elas- 
ticity vs. atomistic elasticity descriptions can affect the ensuing electronic structure of the 
quantum dot. Most notably, the real point group symmetry of the square pyramid is C*2, 
but continuum elasticity spuriously produces a higher C 4 symmetry. 

Regarding the quantitative effects, there are different levels of approximation for coupling 
the strain to the electronic structure. The most general and accurate electronic structure 
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approach is atomistic (e.g., pseudopotentials, tight-binding). There, the full set of atomic po- 
sitions affects the electronic structure. In more approximate electronic structure approaches, 
such as continuum effective-mass, only some aspects of the full, position-dependent-strain 
tensor, e^r) is "felt" by the electronic structure. In these approaches, one considers strain- 
modified potential wells as barriers. Since experiments typically measure electronic energies 
rather than strains, it is instructive to examine these effects. 

Assuming decoupled conduction and valence bands the strain-modified confinement po- 
tential of the conduction-band state is 



£ c (r)=£ c °(r)+a c (r)Tr[e(r)]. 



(23) 



Here, E®(r) is the energy of the conduction-band minimum of the bulk material at r and 
a c , the deformation potential of the conduction-band under hydrostatic deformation. The 
"strain" Hamiltonian of the valence states0 is 



H v = a v (r)Tr(e(r))-b 
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(24) 




-1 
0-10 

where a v is the hydrostatic deformation of the valence states and b and d are uniaxial defor- 
mation potentials for (001) strain and (111) strain, respectively. The effective confinement 
potentials of the valence states are obtained by diagonalizing the strain Hamiltonian coupled 
with the spin-orbit HamiltonianS. Along the [001] direction (z— axis) through the pyramidal 
tip, the shear strains (off-diagonal terms of the strain tensor) are zero and e xx = e yy , and 
thus the effective confinement potentials can be simplified as 



E hh = E° v + a v Tr(e) + -[A so + A s (e)] 



E lh = E° v + a v Tr(e) - ^[A so + A s (e)} + iJ(A- + A*(e)) 2 - ^A-A-(e) 



(25) 
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where A so is the spin-orbit splitting and A s = — 3b[e zz — (e xx + e yy )/2]. 

Figure 7 shows the effective confinement potentials of the conduction and valence-band 
states along the z— axis through the tip of the pyramid. The strain profiles obtained by the 
continuum elasticity and atomistic elasticity are used for the calculation with the same ma- 
terial parameters given in Table ifi Again, the largest difference in confinement potentials is 
found at the interfaces at about 100 meV for the conduction band and 200meV for the valence 
band. The average difference of the confinement potentials inside the dot is about 20 meV 
for the conduction-band state. Although the differences in the strain-modified-confinement 
potentials are small, the band edge states are expected to show different characteristics 
depending upon which strain profile is used for the electronic structure calculation. 

IV. CONCLUSIONS 

We compare the strain distribution of the pyramidal InAs dot grown on a GaAs substrate 
calculated using continuum elasticity and atomistic elasticity. We find a significant difference 
in the strain around the dot interfaces and inside the dot, while the difference in the barrier 
(GaAs substrate and capping layer) is very small. The difference between the two results 
is attributed to the large strain outside the linearity regime of CE, and to the loss of the 
correct atomic symmetry by the CE. 

Acknowledgements This work was supported by United States Department of Energy 
- Basic Energy Sciences, Division of Materials Science under contract No. DE-AC36- 
83CH10093. 
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TABLES 



TABLE I. Ideal bond lengths (d°), elastic constants, and force constants (a and (3) of bulk 
GaAs and InAs. Elastic constants of valence force field method are evaluated by Eq. (8) using a 
and (5 given below. 





4 


C u 


Cl2 


C44 


a 


P 


Material 


(A) 




10 n dyne/cm -1 




10 3 dyne 




GaAs (valence force field) 


2.448 


12.03 


5.70 


5.20 


41.49 


8.94 


GaAs (experimental) 


2.448 


12.11 


5.48 


6.04 






InAs (valence force field) 


2.622 


8.53 


4.90 


3.14 


35.18 


5.49 


InAs (experimental) 


2.622 


8.329 


4.526 


3.959 
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TABLE II. Material parameters used for Fig. 7. All the numbers are given in eV from Ref. [2]. 





A so 


V° 

" V 


a v 


b 


V° 

y c 


a c 


GaAs 


0.34 





1.16 


-1.6 


1.52 


-8.33 


InAs 


0.38 


0.25 


1.00 


-1.8 


0.66 


-6.08 
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FIGURES 

FIG. 1. Schematic diagram of a square pyramidal InAs dot on (001) GaAs substrate. The 
wetting layer (WL) consists of 1 monolayer (ML) of In atoms. Three principal directions, [100], [010] 
and [001] are denoted as X, Y, and Z. The orientation of the pyramidal base isIxY and the 
ratio of the base length (6) and the height (h) is 2 with {110} (grey) and {110} facets. Although 
not shown in the figure, the pyramidal InAs dot is capped by GaAs. 

FIG. 2. Schematics to illustrate how the local strain is calculated. For a cation (Ga or In), three 
vectors ({R}) forming a distorted tetrahedron after atomic relaxation are related to the equivalent 
vectors ({R }) of an ideal tetrahedron via the strain tensor. 

FIG. 3. Strain profiles along the [110] direction at z = h/3 from the base of the pyramid. The 
solid lines are the strain profiles obtained by atomistic elasticity and dotted lines those by contin- 
uum elasticity. The positive and negative signs of X— axes denote the [110] and [110] direction, 
respectively. 

FIG. 4. Solid lines denote Aeij = e^-(CE) — ejj(AE), the difference of each strain component 
obtained by the continuum elasticity and the atomistic elasticity calculations. In III-V Zincblend 
semiconductors, the [110] and [110] directions are inequivalent and, therefore, the symmetry of the 
pyramid is Ci- The C2 symmetry is apparently seen by the difference of the strains (dashed lines) 
along the [110] and [110] directions (see Fig. 1) 

FIG. 5. Strain profiles and the differences along the Z direction through the pyramidal tip. 
The differences between the CE and AE are given on the right-hand side. The discrepancy is 
largest around the interfaces, while the strains in the barrier (GaAs substrate and capping layer) 
agree well within 0.5%. A significant difference is also found inside the quantum dot where the 
InAs experience large compressive strains at about 7 % due to the lattice mismatch. 
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FIG. 6. Relaxation of c eq of GaAs under a biaxial strain obtained by atomistic elasticity is 
compared to the prediction of the continuum elasticity: (a) G = [001] and ej_/eii = — 2C12/C11; 
and (b) G = [110] and e ± /e\\ = -(C n + 3C 12 - 2C AA )/{Cu + C 12 + 2C 44 ). At the infinitely 
small limit, the —e±/e\\ of the atomistic elasticity (solid lines) coincides with that of the continuum 
elasticity (dashed lines). The discrepancy between the CE and AE increases for the larger biaxial 
strains. 

FIG. 7. (a) Confinement potentials by Eq. (15) and (16), along Z direction through the pyra- 
midal tip in Fig. 1, predicted by the strain profiles obtained by the atomistic elasticity. All the 
energies are measured with respect to the valence band maximum of the bulk GaAs at equilibrium, 
(b) The difference of the confinement potentials by continuum elasticity and atomistic elasticity. 
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